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Abstract 

Inclusive K^K^ production in deep inelastic ep scattering at HERA has been 
studied with the ZEUS detector using an integrated luminosity of 120 pb~^. 
Two states are observed at masses of ISSTlg MeV and 1726 ± 7 MeV, as well as 
an enhancement around 1300 MeV. The state at 1537 MeV is consistent with the 
well established /2(1525). The state at 1726 MeV may be the glueball candidate 
/o(1710). However, it's width of 38^14 MeV is narrower than 125 ± 10 MeV 
observed by previous experiments for the /o(1710). 
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1 Introduction 



The K^K^ system is expected to couple to scalar and tensor glueballs. This has motivated 
intense experimental and theoretical study during the past few years |2I21- Lattice QCD 
calculations jH] predict the existence of a scalar glueball with a mass of 1730 ±100 MeV 
and a tensor glueball at 2400 ± 120 MeV. The scalar glueball can mix with qq states with 

1 = from the scalar meson nonet, leading to three J^*" = O"^"*" states, whereas only two 
can fit into the nonet. Experimentally, four states with J^*^ = O"*"^ and / = have been 
established g]: /o(980), /o(1370), /o(1500) and /o(1710). 

The state most frequently considered to be a glueball candidate is /o(1710) but its 
gluon content has not yet been established. This state was first observed in radiative J/ ijj 
decays j^I and its angular momentum J = was established by the WA102 experiment 
using a partial- wave analysis in the K~^K~ and K^K^ final states |n|. A recent publication 
from L3 [2j reports the observation of two states in 77 collisions above 1500 MeV, the 
well established /2(1525) j3] and a broad resonance at 1760 MeV. It is not clear if this 
state is the /o(1710). Observation of /o(1710) in 77 collisions would indicate a large quark 
content. 

The ep collisions at HERA provide an opportunity to study resonance production in a 
new environment. The production of K'^ has been studied previously at HERA jHlini- In 
this paper, the first observation of resonances in the K^K^ final state in inclusive ep deep 
inelastic scattering (DIS) is reported. 

2 Experimental set-up 

The data were collected by the ZEUS detector at HERA during the 96-00 running peri- 
ods. In 96-97, HERA colhded 27.5 GeV positrons with 820 GeV protons. In 98-00, the 
proton energy was 920 GeV and both positrons and electrons were collided with protons. 
The measurements for e^p {e~p) interactions^ are based on an integrated luminosity of 
104 pb"^ (17pb-i). 

A detailed description of the ZEUS detector can be found elsewhere ^Uj. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) ^T], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 

^ Hereafter, both e+ and e~ arc referred to as electrons, unless explicitly stated otherwise. 
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polar-angle^ region 15° < 9 < 164°. The relative transverse- momentum resolution for full- 
length tracks is a{pT)/pT = 0.0058pT© 0.0065 ©0.0014/pt, with pr in GeV. The tracking 
system was used to establish the primary and secondary vertices. 

The high-resolution uranium-scintillator calorimeter (CAL) consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/v^ for electrons 
and a{E)/E = 0.35//E for hadrons, with E in GeV. 

The energy of the scattered electron was corrected for energy loss in the material be- 
tween the interaction point and the calorimeter using a small-angle rear tracking detector 
(SRTD) [131111 and a presampler (PRES) [13111. 

3 Kinematic reconstruction and event selection 

The inclusive neutral current DIS process e{k) + p{P) e(/c') + X can be described 
in terms of the following variables: the negative of the invariant-mass squared of the 
exchanged virtual photon, = —q^ = — [k — k'Y\ the fraction of the lepton energy 
transferred to the proton in the proton rest frame, y = {q ■ P)/ {k ■ P); and the Bjorken 
scaling variable, x = Q"^ /{2P • q). 

A three-level trigger system was used to select events online jTUj. The inclusive DIS 
selection was defined by requiring an electron found in the CAL. In certain run periods, 
corresponding to 83% of the total luminosity, the inclusive selection was not available 
for low {Q"^ < 20 GeV^) events. For these periods, an additional selection with a 
requirement of at least one forward jet identified with the kt algorithm and having 
transverse energy Et > 3 GeV {E^ > 4 GeV in the 2000 running period), pseudorapidity 
< ?7jct < 3 (1.5 < ?7jct < 3.5 in the 1996-1997 running period) was used. 

The DIS offline event selection was based on the following requirements: 

• a primary vertex position, determined from the tracks fitted to the vertex, in the range 
I .^vertex I < 50 cm, to rcducc the background events from non-ep interactions; 

• Ee > 8.5 GeV, where Ee is the energy of the scattered electron reconstructed in the 
calorimeter; 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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• 42 < 6 < 60 GeV, where 6 = — cos6'j), Ei is the energy of the i^^ 
calorimeter cell and 6i is its polar angle as viewed from the primary vertex. The sum 
runs over all cells. This cut further reduces the background from photoproduction and 
events with large QED initial-state radiation; 

• Uc ^ 0.95, to remove events with misidentified scattered DIS electrons; j/e is the value 
of y reconstructed using the scattered electron measurements; 

• yjB > 0.04, to remove events with low hadronic activities; ?/jb is the value of y 
reconstructed using the Jaquet-Blondel method [T7j : 

• the position of the scattered lepton candidate in the RCAL was required to be outside 
a box of ±14 cm in X and Y, which corresponds approximately to 6*010 ~ 176°, where 
6'eie is the polar angle of the scattered electron; 

• a maximum of 40 tracks per event. This cut reduces the background from false 
pair candidates, removing only 7% of the DIS events. 

The non-ep and photoproduction background in the selected sample was negligible. The 
data were not corrected for the biases introduced by the trigger requirements and selection 
cuts. 



4 Selection of i^^-pair candidates 

The Kg meson candidates were reconstructed using tracks with at least 38 hits in the 
CTD and pseudorapidity within ±1.75. The pseudorapidity cut and the required min- 
imum number of hits ensured a good momentum resolution and a minimum transverse 
momentum of 0.1 GeV. In each event, oppositely charged track pairs assigned to a sec- 
ondary vertex were combined to form candidates. Both tracks were assigned the mass 
of a charged pion and the invariant-mass M(7r^7r~) calculated. The secondary vertex 
resolution for these events, estimated using MC studies, is 2 mm in X and Y, and 4 mm 
in Z. 

Additional requirements were applied to the selected candidates: 

• Pt{Ks) > 200 MeV, for each K° candidate; 

• 2 < d < 30 cm, where d is the decay length of the candidate; 

• dxY < 4 mm and dz < 5.5 mm, where dxY and dz are, respectively, the projections 
on the XY plane and Z axis of the vector defined by the primary interaction point and 
the point of closest approach of the K'^ candidate; 
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• (^XY < 0.12, where 6xy is the (colhnearity) angle between the candidate momen- 
tum vector and the vector defined by the interaction point and the decay vertex 
in the XY plane; 

• > 110 MeV, for each candidate, where the Armenteros-Podolanski variable 

is the projection of the candidate pion momentum onto a plane perpendicular to the 
candidate line of flight [TH] . 

The cuts on the decay length, distance of closest approach and colhnearity angle signifl- 
cantly reduce the non-Kg background as determined by Monte Carlo (MC) simulations. 
After the cut, backgrounds from A, A and photon conversions are negligible. Only 
events with at least two selected candidates were kept for further analysis. 

Figure Q shows the invariant-mass distribution for candidates in the range 0.45 < 
M{7i~^7r~) < 0.55 GeV after the Kg-paii candidate selection. The distribution was fltted 
using one linear and two Gaussian functions. The linear function flts the background, one 
of the Gaussians flts the peak region in the central 7T~^7r~ invariant mass distribution, and 
the other Gaussian improves the flt at the tails. The two Gaussians were constrained to 
have the same mean value. A mass of 498 MeV and standard deviation width of 4 MeV 
were obtained from the flt with the central Gaussian, and a standard deviation width of 
8.5 MeV was obtained from the flt with the other Gaussian. The normalization factor 
between the narrower and broader Gaussians is approximately 5.5. The invariant-mass 
width is dominated by the momentum resolution of tracks reconstructed with the CTD. 
Only Kg candidates in the region of ±10 MeV around the fltted central mass were used 
to reconstruct the K^K^ invariant-mass. High-statistics samples of Monte Carlo events, 
generated without resonances, were used to conflrm that the event-selection criteria did 
not produce artiflcial peaks in the K^K^ invariant-mass spectrum. 

Figure El shows the distribution in x and of selected events containing at least one 
pair of candidates. The kinematic variables were reconstructed using the double 
angle method ^H]- The virtual-photon proton centre-of-mass energy was in the range 
50<W < 250 GeV. 

5 Results 

The K^Kg spectrum may have a strong enhancement near the K^K^ threshold due to 
the /o(980)/ao(980) state [213121] • Since the high K^K^ mass is the region of interest for 
this analysis, the comphcation due to the threshold region is avoided by imposing the cut 
cosOkoro < 0.92, where Okoro is the opening angle between the two candidates in the 
laboratory frame. 
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After applying all selections, 2553 -R'^-pair candidates were found in the range 0.995 < 
M{K^^K^^) < 2.795 GeV, where M{K^K^) was calculated using the mass of 497.672 MeV 

The momentum resolution of the CTD leads to an average M{K^K^) resolution which 
ranges from 7 MeV in the 1300 MeV mass region to 10 MeV in the 1700 MeV region. Fig- 
ure El shows the measured K^K^ invariant-mass spectrum. Two clear peaks are seen, one 
around 1500 MeV and the other around 1700 MeV, along with an enhancement around 
1300 MeV. The data for cos^i^oxo > 0.92 are also shown. The mass scale uncertainties in 
the region of interest, arising from uncertainties in the magnetic field, is at the per mille 
level. 

The distribution of Fig.Olwas fitted using three modified relativistic Breit-Wigner (MRBW) 
distributions and a background function U{M); 

3 p 

F(M) = Y(—, "^to 9 9 ) + U(M) , (1) 

where j is the effective resonance width, which takes into account spin and large width 
effects [22]; m^^i is the resonance mass, and M is K^K^ invariant mass. The background 
function is 



U{M) = A-{M- 2mKof ■ e^^V^^™^ , (2) 

where A, B and C are free parameters and nixo is the mass defined by Hagiwara et 
al j3] . Monte Carlo studies showed that effects of the track-momentum resolution on the 
mass reconstruction were small compared to the measured widths of the states. Therefore 
the resolution effects were ignored in the fit. 

Below 1500 MeV, a region strongly affected by the cosOkok^ cut, a peak is seen around 
1300 MeV where a contribution from /2(1270)/a2(1320) is expected. This mass region 
was fitted with a single Breit-Wigner. 

Above 1500 MeV, the lower-mass state has a fitted mass of 1537t^ MeV and a width 
of b^tlt MeV, in good agreement with the well established /2(1525). The higher-mass 
state has a fitted mass of 1726 ± 7 MeV and a width of SSl^" MeV. The widths reported 
here were stable, within statistical errors, to a wide variation of fitting methods including 
those using 30 MeV bins rather than the default 15 MeV bins. The width is narrower 
than the PDG value (125 ± 10 MeV) H] reported for /o(1710), but when it is fixed to this 
value, the fit is still acceptable. For the purposes of the following discussion, this state is 
referred to as /o(1710). 

The masses, widths and number of events from the fit with statistical errors are given 
in the top row of Table H It should be noted that there are correlations between the 
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obtained parameters. In particular the width and the number of events for each state are 
highly correlated. The mass spectrum is also consistent with the background function at 
higher masses, but masses above 2795 MeV were not included in the fit due to limited 
statistics. The sensitivity of the data to the widths of the resonances was studied. Several 
fits were performed fixing the width of the states /2(1525) and /o(1710) to their PDG 
values; the results are shown in rows 2 to 4 of Table 

In the literature, there are several states reported in the mass region near 2000 MeV 
m, namely /2(1950) and /o(2020) which need confirmation, and /2(2010), 04(2040) and 
/4(2050) which have been confirmed. While no visible structure in this mass region exist 
in the data, these states may affect the background shape. The effect of the inclusion 
of a state in this region was examined, with the result shown in row 5 of Table H The 
description of the data between 1800 and 2000 MeV is improved with respect to the other 
fits. However, there is no improvement in the overall /dof . 

The K^K^ spectrum after all cuts was also fitted with the background function only. The 
fit can be rejected with a 99.4% confidence level using a simple test over the mass 
region 1000 to 2800 MeV. 

It was found that 93% of the i^°-pair candidates selected within the detector and trigger 
acceptance are in the target region of the Breit frame j2Hj, the hemisphere containing the 
proton remnant. Of the /C^-pair candidates in the target region, 78% are in the region 
Xp = 2pb/Q > 1, where pb is the absolute momentum of the K^K^ in the Breit frame. 
High Xp corresponds to production of the /T^-pair in a region where sizeable initial state 
gluon radiation may be expected. This is in contrast to the situation at e^e~ colliders 
where the particles entering the hard scattering are colourless. 



6 Conclusions 



The first observation in ep deep inelastic scattering of a state at 1537 MeV, consistent 
with /2(1525), and another at 1726 MeV, close to /o(1710), is reported. There is also 
an enhancement near 1300 MeV which may arise from the production of /2(1270) and/or 
02(1320) states. The width of the state at 1537 MeV is consistent with the PDG value for 
the /2(1525). The state at 1726 MeV has a mass consistent with the glueball candidate 
/o(1710), and is found in a gluon-rich region of phase space. However, it's width of 381^° 
MeV is narrower than the PDG value of 125±10 MeV for the /o(1710). 
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Fit 




/2(1270)/a^(1320) 


/2(1525) 


/o(1710) 


/j(1980) 






mass 


width 


events 


mass 


width 


events 


mass 


width 


events 


mass 


width 


events 


1 


0.97 




244^1 


414lli 


15371^ 


501^^ 


841^} 


1726 ± 7 


38l?o 


7411 








2 


0.96 


1272 ± 16 


24011 


42011^7 


1539 ± 10 


76 


107 ± 30 


1727 ± 7 


39 ±20 


7(^+28 








3 


1.02 


1276 ± 16 


2581^° 


48011^0 


1536 ± 8 


49i? 


oc+38 
»t)_27 


1726 ± 13 


125 


122 ± 40 








4 


1.02 


1274 ± 15 


251^1 


476111 


1538 ± 10 


76 


108li 


1728 ± 13 


125 


1201^^ 








5 


1.00 


1283 ± 15 


26011 


5061218 


154011^ 


7Q+43 


116li^ 


1727 ± 7 


47+23 


911^^ 


19701^^ 


138lL^^ 


74 ±40 


Particle Data Group 2002 Values (MeV) [H 






1275.5 ± 1.2, 185l:|^ 
1318 ±0.6, 104.7 ± 1.9 


1525 ± 5, 76 ± 10 


1713 ± 6, 125 ± 10 





Table 1: Normalized , masses in MeV, widths in MeV, and number of events 
extracted from the K^K^ invariant mass fits (errors are statistical only.) Widths 
reported without errors were fixed to the PDG values listed in the last row. 
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Figure 1: The distribution of tt^tt^ invariant mass for events with at least two 
candidates passing all selection cuts. The solid line shows the result of a fit 
using one linear and two Gaussian functions. 
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ZEUS 




X 



Figure 2: The distribution in x and of events passing all selection cuts. The 
dashed lines delineate approximately the kinematic region selected. The solid line 
indicates the kinematic limit for HERA running with 920 Ge V protons. 



11 



ZEUS 




M(K»K°) (GeV) 



Figure 3: The K^K^ invariant-mass spectrum for pair candidates with 
cosOkoko < 0.92 (filled circles). The thick solid line is the result of a fit using three 
Breit-Wigners (thin solid lines) and a background function (dotted- dashed line). 
The Kg pair candidates that fail the cosOkoko < 0.92 cut are also shown (open 
circles). 
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